Ab initio studies of the tunneling magneto-Seebeck effect: influence of magnetic 
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We found a strong influence of the composition of the magnetic material on the temperature 
dependence of the tunneling magneto-Seebeck effect in MgO based tunnel junctions. We use ab 
initio alloy theory to consider different Fe x Co\- x alloys for the ferromagnetic layer. Even a small 
change of the composition leads to strong changes in the magnitude or even in the sign of the 
tunneling magneto-Seebeck effect. This can explain differences between recent experimental results. 
In addition, changing the barrier thickness from six to ten monolayers of MgO leads also to a non- 
trivial change of the temperature dependence. Our results emphasize that the tunneling magneto- 
Seebeck effect depends very crucially and is very sensitive to material parameters and show that 
further experimental and theoretical investigations are necessary. 
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The recently theoretically predicted [l| and experimen- 
tally confirmed @, 0] tunneling magneto-Seebeck (TMS) 
effect in MgO based tunnel junctions belongs to the 
new field of spin caloritronics 0, IH- In this field the 
spin-dependent charge transport is combined with en- 
ergy or heat transport. This means that the spin de- 
gree of freedom is exploited in thermoelectrics [6]. Be- 
sides the TMS currently investigated effects in the field 
of spin caloritronics are the spin-Seebeck effect 0, 3, 
the magneto-Seebeck effect in metallic multilayers 9j, 
the thermal spin-transfer torque [Icj , the spin-dependent 
Seebeck effect [ll[, thermally excited spin-currents 12 1, 
and the magneto-Peltier cooling [13]. In this letter we 
investigate the role of the ferromagnetic lead material on 
the TMS. We show that not only the temperature de- 
pendence but even the sign of the TMS effect depends 
crucially on the material composition. 

The TMS effect is the change of the Seebeck coefficient 
with a change of the magnetic orientation of the ferro- 
magnetic leads relative to each other in a tunnel junction. 
The size of the effect is given by the TMS ratio 
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where S p (S AP ) is the Seebeck coefficient for parallel 
(anti-parallel) magnetic orientation of the ferromagnetic 
leads. Therefore, the TMS effect is similar to the tunnel- 
ing magnetoresistance effect (TMR) [l4, 15 1, where one 
considers the change of the electrical resistance with a 
change of the magnetic orientation. Note, that in con- 
trast to the resistance the Seebeck coefficient can be pos- 
itive or negative. Hence, the TMS ratio can have diver- 
gences whenever one of the Seebeck coefficients in Eq. ([T]) 
is zero. 

The tunnel junctions we study consist of a 
Fe x Coi- x (Q0l)/MgO/Fe x Coi -^(OOl) structure embed- 
ded between semi-infinite leads. These leads are just act- 
ing as reservoirs and are modeled by Cu in the hcc-Fe 



structure. For the Fe x Co\- x alloy we use a fixed lattice 
constant of 0.287nm for all compositions. The thickness 
of both ferromagnetic leads is 20 monolayers and only 
symmetric junctions are considered. For the MgO bar- 
rier we use six and ten monolayers. Like in our previous 
studies [l|, [H we use the ideal positions for the inter- 
face layers. In particular, no relaxation effects at the 
Fe x Co\- x jMgO interface are taken into account. That 
way, we really focus on the influence of the change of 
the electronic structure in the lead material by alloying. 
Although in most cases, e.g. by sputtering techniques, 
FepCoi- x is only stable in bcc structure for x > 0.3 
j 1 €31] we consider the whole concentration range. Actu- 
ally, with methods like molecular beam epitaxy it is also 
possible to grow pure Co leads in bcc-structure [17 j . 

For the description of the Fe x Co\- x alloy we em- 
ploy the coherent potential approximation (CPA) [lj| [l9[ 
recently implemented in our Korringa-Kohn-Rostoker 
(KKR) method Within the CPA the alloy is de- 

scribed by an effective medium, which is calculated self- 
consistently. For the description of the transport prop- 
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erties so-called vertex corrections are essential 
The CPA together with the vertex corrections leads ba- 
sically to the same result as the supercell approach, where 
one has to average the transport properties over a larger 
number of different supercells [20j . In comparison to the 
supercell approach, the advantage of the CPA is a lower 
computational effort and the possibility to use an arbi- 
trary composition. The transport coefficients, in particu- 
lar the transmission function T(E), are calculated using 
the non-equilibrium Green's function formalism imple- 
mented in the KKR method [3] including vertex cor- 
rections for the CPA [20j]. From T{E) we calculate the 
moments 
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where f(E, n, 0) is the Fermi occupation function at a 
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given energy E, electrochemical potential fi, and temper- 
ature O. In linear response the conductance G and the 
Seebeck coefficient S are given by [25[ 
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These quantities are calculated for parallel and anti- 
parallel magnetic orientation of the Fe x Co\- x layers to 
eventually calculate the TMS ratio according to Eq. (pj . 
Note, that the conductance is basically the area under the 
transmission function times the derivative of the occupa- 
tion function, whereas the Seebeck coefficient is propor- 
tional to the expected value (or first moment) of the very 
same product. Therefore, the Seebeck coefficient is de- 
termined by the asymmetry of the transmission function 
with respect to the Fermi level. Normally, the geometric 
mean of a function is more affected by small modulations 
of the function in comparison to the area under that func- 
tion. Consequently, it can be expected that the Seebeck 
coefficient can depend crucially on changes of the trans- 
mission function e.g. due to alloying of the magnetic 
material. All calculations are done in the atomic sphere 
approximation and the cut-off for the angular momen- 
tum is 3. The k-point grids consist of 576 and 40,000 
points in the whole Brillouin zone for the self-consistent 
and transport calculation, respectively. 

One of our main results is shown in Fig. [T] There, the 
TMS ratio at room temperature for two different barrier 
thicknesses is given as a function of the Fe x Co\- x com- 
position. It clearly shows that not only the magnitude 
but also the sign of the TMS ratio is very sensitive to the 
actual composition. This means that even small changes 
in the composition can drastically change the magnitude 
and even the sign of the TMS ratio. This could explain 
the small experimentally observed values 0, Q in com- 
parison to the large values predicted for the pure materi- 
als [l| , because in experiments usually compositions close 
to fifty-fifty Fe and Co are used. In addition, differences 
in the composition of the magnetic material can be a rea- 
son for the different signs observed in the experiments in 
Ref . | and Ref . & 

Increasing the barrier thickness from six to ten mono- 
layers can change the values significantly depending on 
the composition. In particular, the TMS ratio can in- 
crease as well as decrease with increasing barrier thick- 
ness and even the sign can change. To analyze the TMS 
further the bottom viewgraph of Fig. [1] shows the See- 
beck coefficients at room temperature for parallel and 
anti-parallel magnetic configuration that leads to the dis- 
cussed TMS ratios according to Eq. (JTJ) . Likewise, here 
no general trend is visible and the Seebeck coefficient can 
be decreased or increased by increasing the barrier thick- 
ness. 

Up to now we discussed the TMS only at room tem- 
perature. In Fig. [5] we show for the different considered 
compositions the temperature dependence of the TMS 
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FIG. 1: (Color online) Top: Tunneling magneto-Seebeck 
(TMS) ratio as a function of Fe x Coi- x composition at room 
temperature for two different barrier thicknesses. Bottom: 
Seebeck coefficient for parallel S p and anti-parallel S AP mag- 
netic configuration of the ferromagnetic layers as a function 
of the Fe x Coi- x composition at room temperature for two 
different barrier thicknesses. 



ratio for six and ten monolayers of MgO. Note, that 
the temperature enters the calculation only in the oc- 
cupation function in Eq. @. The temperature depen- 
dencies for the pure materials were already discussed in 
Ref. [l|. Fig. [5] shows that the temperature dependence 
is rather involved and that there is also no clear trend 
visible when the composition is changed. We illustrate 
this with two examples. First, in the lower left panel all 
shown compositions have a moderate increase except for 
the Fcq.^Coq^ alloy, which has a strong increasing TMS 
ratio with temperature. Second, in the lower right panel 
for FeojCoo.3 and Feo.gCon.i the TMS ratio is decreas- 
ing with increasing temperature but for the composition 
in between Feo.sCoo.2 the TMS ratio is increasing. Be- 
sides the dependence on the composition also the change 
in the barrier thickness can lead to quite different temper- 
ature dependencies, for example visible for -Fen^Con.s, 
Fe . 5 Co Q . 5 , and Fe . s Co .2. 

Particularly striking features are divergences of the 
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FIG. 2: (Color online) Tunneling magneto-Seebeck (TMS) ratio as a function of temperature for pure Fe and Co (as already 
published in Ref. [l|) and for different alloy compositions of Fe x Co\- x . For the barrier we use six and ten monolayers of MgO. 



TMS ratio which occur when one of the Seebeck coef- 
ficients goes through zero. Besides for pure Fe only for 
the -Feo.iCorj.g alloy a divergence is observed. Another 
interesting feature is a sign change of the TMS ratio with 
temperature, which is also observed experimentally 
To get an overview of when the sign changes occur we 
compile in Tab. U the temperatures at which the sign 
changes occur and in which direction it proceeds. For 
the compositions where no sign change occurs we give 
the corresponding sign of the TMS ratio over the whole 
temperature range. We see that there is a sign change of 
the TMS ratio only for certain compositions and that it 
also depends on the barrier thickness. Moreover, a com- 
parison of Feo.eCooA and Feo.yCoo.s shows again that 
the sign of the TMS ratio crucially depends on the com- 
position of the ferromagnetic material. 

For completeness, we give in Fig. [3] the corresponding 
temperature dependencies of the Seebeck coefficients for 
parallel and anti-parallel magnetic configuration, which 
lead to the temperature dependencies of the TMS shown 
in Fig. [2] It shows that a vanishing Seebeck coefficient 
in the anti-parallel configuration is responsible for all the 
observed divergencies, which is just by accident and has 
no physical reason. 



Up to now, all the presented results seem unsystem- 
atic with respect to the alloy composition. The underly- 
ing physical quantity is the transmission function T(E), 
which determines the transport properties according to 
Eqs. ([2|) and ([3]). Consequently in Fig. |4] we look at 
T(E) for the parallel and anti-parallel configuration for 
the two considered barrier thicknesses. For the transmis- 
sion in the parallel configuration (left panels in Fig. @| 
a systematic and continuous change with the ferromag- 
netic layer composition is visible. For example, for six 
monolayers of MgO there is a peak in T P (E) below the 
Fermi level, which moves continuously to higher energies 
with increasing Fe concentration. For an Fe concentra- 
tion of about 70% a second peak occurs at energies below 
the Fermi level. One has to keep in mind that the See- 
beck coefficient is proportional to the first moment of the 
transmission function times the derivative of the occupa- 
tion function (see Eqs. ^ and |(3J|). Consequently, only 
the asymmetry of T(E) is important and this is changed 
by the discussed movement of peaks. With this the given 
temperature dependencies of S in Fig. [3] can be under- 
stood. Moreover, due to the continuous change of T P (E) 
one can interpolate to other compositions without doing 
a full ab initio calculation. 
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TABLE I: Temperatures of the sign changes T ac of the tunneling magneto-Seebeck (TMS) ratio, if present, extracted from the 
different temperature dependencies given in Fig. [2] In addition, the direction of the sign change is given in the table: n — > p 
(p — ► n) means a change from negative to positive (positive to negative) with increasing temperature. If no sign change is 
present n or p gives the sign of the TMS ratio over the whole temperature range. 
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FIG. 3: (Color online) Seebeck coefficients for parallel (black) and anti-parallel (red) magnetic configuration for different 
Fe x Co\- x compositions. For the barrier we use six (solid) and ten (dashed) monolayers of MgO. 



The reason for this continuous change in the parallel 
configuration can be explained with the band structure. 
First, it is important to realize that the transport is dom- 
inated by the majority spin in Fe and Co where a Ai 
band is present and dominates the transport [26l - l28| . By 
alloying this band is only weakly affected and the com- 
position primarily shifts the position of the Fermi level. 
And this shift of the Fermi level corresponds to the shift 
of the transmission function with the composition. 

However, this is no longer the case for the anti-parallel 
configuration. Here, the simple Ai band does not con- 
tribute to the transport and the transport is dominated 
by pockets within the Brillouin zone [26J. This is clearly 
visible in the rich structure of T AP for the pure mate- 
rials. In the alloy case the band structure is broadened 
by disorder-scattering leading to a "washed out" trans- 
mission function. Going through the different composi- 
tions no clear continuous change is visible, in particular, 
close to the pure materials. Therefore, the rather in- 
volved temperature dependencies of the TMS ratio can 
be traced back to the complicated change with compo- 
sition of the transmission function in the anti-parallel 
magnetic configuration. 

In conclusion, we show that the TMS ratio is cru- 
cially dependent on the alloy composition of the magnetic 
material. The behavior seems unsystematic but can be 
traced back to a simple concentration dependence of the 



parallel transmission and a complicated change of the 
transmission function in the anti-parallel magnetic con- 
figuration. This leads to a strong dependence not only of 
the magnitude but also of the sign of the TMS ratio even 
for small changes in compositions. In general, the TMS 
ratios are smaller for compositions close to fifty-fifty in 
comparison to compositions closer to the pure materials. 
Moreover, small changes in compositions can cause the 
different signs of the TMS, which were observed in exper- 
iments 0, y| . To achieve a full understanding of the TMS 
in the MgO based tunnel junctions further experimental 
and theoretical studies are necessary. 
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